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By considering full-field string solutions of the Abelian-Higgs model, we modify the model of a fluid 
of strings (which is composed of Nambu strings) to obtain a model for a "fluid of vortices." With this 
model, and following closely Soleng's proposal of a fluid of strings as the source of a Milgrom-type 
correction to the Newton dynamics, we determine quantitatively the modified dynamics generated 
by a static, spherical fluid of vortices. 
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Since its development more than eighty years ago, Ein- 
stein's theory of General Relativity has enjoyed many 
successes in describing gravity at very different scales, 
from the orbit of Mercury to the prediction of black holes 
or the Big Bang model. The theory, however, encoun- 
ters serious problems when confronted to the motion of 
galaxies and galaxy clusters, where it seems to imply that 
there is more mass than is observed. This is the so-called 
"missing mass" problem, to which several solutions are 
being investigated: first, General Relativity may have to 
be modified at such distances; to conform to the obser- 
vations, the modified theory would then have to yield a 
gravitational acceleration decreasing as 1/r. Second, the 
Universe may be filled with dark matter, which would 
act gravitationally but would not be observable. Third, 
non-gravitational forces may play an important role at 
very large scales. Finally, of course, it might also prove 
necessary to invoke a combination of several of these so- 
lutions. 

There are good reasons, however, to believe that Gen- 
eral Relativity is not the definitive theory of gravitation, 
in that it is not a quantum theory. Currently, string/M- 
theory is seen by many as the best candidate to unify 
the four fundamental forces, and it is well known that 
in such theories Einsteinian gravity must be replaced by 
more complicated scalar-tensor models. In low-energy 
M-theory, for instance, one must compactify seven of the 
eleven spacetime dimensions to obtain General Relativ- 
ity from eleven-dimensional supergravity. Although it is 
widely believed that the extra dimensions only play a 
significant role at very short distances, there exist mod- 
els where gravity becomes higher dimensional jjj at very 
large distances as well. 

A rather simple, but effective, approach to the prob- 
lem is to modify the Newton force at large distance. Mil- 
grom H proposed to write the real gravitational acceler- 
ation g as a function of the Newtonian acceleration gN, 



/' I — ) g = g.v 



(1) 



where ao is a constant, and n(x) ~ 1 for x <C 1 and 
p(x) x for i>1. This modified Newtonian dynamics 
was used to explain, without the need for any dark mat- 
ter, the observed gravitational behavior of galaxies and 



galaxy clusters P,f|]. In particular, it yields the a con- 
stant velocity curve at large galactic radius, V 4 ~ MGao, 
which implies that the true acceleration at large distances 
is proportional to 1 /r rather than 1 jr 2 . 

The main problem with this model is that it is sup- 
ported only by its phcnomcnological success and has no 
theoretical basis. 

By using a model of a fluid of ordered strings, devel- 
oped some time ago by Letelier [||6), and usual General 
Relativity, Soleng found that a perfect fluid surround- 
ing a point mass M leads to a force 
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where I is an integration constant. The parameter a 
characterizes the fluid's equation of state, and remains 
unspecified. 

In this Letter, our main goal is to introduce a new 
model describing a "perfect fluid of vortices," which is 
a relative of the well-known string fluid 0, except that 
we replace the Nambu strings by (thick) solutions of the 
full-field equations of the Abelian-Higgs model, namely 
Nielsen-Olesen (NO) vortices ||. Because we choose a 
particle model to describe the vortices, we are able to 
choose more realistically the parameter a appearing in 
Soleng's solution, and therefore also to determine the 
modified dynamics of the perfect vortex fluid surrounding 
a mass M. 

Let us begin by describing the original string fluid 
model introduced in || (see also (|]). The starting point 
is the observation that the natural generalization of the 
perfect fluid energy-momentum tensor to the case of one 
dimension extended objects with equations of state "ten- 
sion equals energy density" is 



(p+p) (u^-x'xn-pg 



(3) 



where p is the string energy-density that is equal to its 
tension, p > (p < 0) is a transversal pressure (tension), 
and u M and are two smooth vector fields satisfying 



u^Xn = 0. 



(4a) 
(4b) 
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At any spacetime point represents the string velocity 
and x* 1 the string direction at each point of the string. 

The metric for a static spherically symmetric space- 
time is 

ds 2 = e 2 ^dt 2 ~ e 2X ^dr 2 - r 2 {d9 2 + sin 2 9d V 2 ) . (5) 

The general solution of the Einstein equations for the 
energy-momentum tensor (^) and metric (|J) was studied 
in raj. For the equation of state 



and 



p = -ap, 
it was found (for a ^ 2) that 
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where £ is an integration constant. This metric repre- 
sents a bunch of radially directed Nambu strings inter- 
acting transversally, in a way described by the equation 
of state ([j|) . 

Our aim now is to model this transversal interaction 
by considering a particle model of the strings in the fluid. 
More precisely, we consider an Abelian-Higgs model with 
matter Lagrangian 

Cm = (V^) (V^) f - \f^ v F^ - ± ($t$ _ jf) 2 , ( 8 ) 

where $ is a complex Higgs field, T> ^ is a gauge-covariant 
derivative defined by = + ieA^Q and F^ v is 

the strength field associated with the gauge potential A,, 
that is, F^ v = 2d [fl A v] . 

We can simplify this theory by making the conven- 
tional Ansatz in cylindrical coordinates {t, z, R, <f\ 



$ = r}X(R)e 



:N<p 



A„ = - (P„ - V aX ) ■ 



(9a) 
(9b) 



We can simplify this further by simultaneously rescal- 
ing our coordinates and P M by the string's width wn = 
l/(y/\rj). This corresponds to measuring distances in 
string rather than Planck units, and can be most easily 
done by letting wjj = 1 everywhere. The only parameter 
left in this theory is the Bogomol'nyi parameter 
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(10) 



The NO solution is found by assuming that N = 
1,P M = P(R)Wfj,(f>, in which case the Lagrangian and 
equations of motion become (after removal of a global 
multiplying constant) 



— Cm — X' 



pi 2 



X 2 P 2 
~R 2 ~ 



Y> Y 2 P 2 1 

X" + — - ±J— - -X (X 2 - 1) = 
R R 2 2 V ; 
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X 2 P = 0, 



(12a) 
(12b) 



where a prime denotes differentiation with respect to R. 
To complement these equations, we must set the appro- 
priate boundary conditions, 



X(oo) = P(0) = 1. (13) 



X(0) = P(oc) = 0, 



With this we can solve the NO equations fll2|) numerically 
for various values of the parameter (3. Sample solutions 
are plotted on figure [l| 
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FIG. 1. The Nielsen-Olesen solution X(R) and P(R) for a 
few values of the Bogomoln'yi parameter, j3 = 1/8, 1/4, 1/2, 
1, 2, 4 and 8. 

The ener gy-m omentum tensor for the NO solution is 
given from (11) by 



T» v = Diag (p,p z ,p fl ,p ) , 



where 



,2 P' 2 

p= x +< 9 #r + 
Pz = -p = A 

p'2 

PR = -x' 2 -p 1 ^ + 

pl2 



X 2 P 2 
~R 2 ~ 

X 2p2 

~rT~ 

X 2 P 2 
~R 2 ~ 



1 , 9 

4 V 



-(X 2 
4 V 

-(X 2 

4 V 



-1 



-1 



(14) 
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Note that for /3 = 1 the equations of motion reduce to 
the Bogomol'nyi equations, 
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(16) 



which yields pr — p^ — 0. For the purpose of finding the 
vortex's equation of state, we must integrate the energy 
and the pressures; these are shown in figure 0. 




P 

FIG. 2. Total energy, radial pressure and azimuthal pres- 
sure (top to bottom) for the NO solution, as functions of f3. 

Finally, we can write the equations of state of the vor- 
tex as 



PR = q-r{P)p 
p e = a g (/3)p. 



(17a) 
(17b) 



The functions an(0) and ag(0) are given on figure^. 

We are now in a position to determine realistic values 
for cur and ag. Because the Abelian-Higgs theory is a 
toy model, there is no experimental determination of the 



coefficients A and e appearing in the definition of f3, but 
we can compute it from the masses of the Higgs and the 
gauge boson. We assume here that 



77 < m H (GeV) < 300 



(18) 



and that the most probable value is tojj ~ 170 GeV [ 
For the gauge boson, we assume a W or a Z°, with 



m w = 80GeV, 
This yields 



m z = 82GeV. 



(19) 



0.88 < < 14.0 (20a) 
-0.011 < a R < 0.239 (20b) 
-1.103 x 10~ 4 < ag < -2.51 x 10~ 5 , (20c) 

the most probable value for tojj implying (3 = 4.4, an — 
0.134 and ag = -4.27 x 10~ 5 . 

We can now apply what we found to a spherical domain 
with metric (^J) containing a fluid of vortices, by writing 
the equation of state for this fluid as 



(21) 



where —5.51 x 10~ 3 < k < 0.119 and most probable value 
K = 0.067. Insertion of this into Soleng's result Q, via 
k = —1/a, gives 



2v -2A -i 

e = e =1 
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(22a) 
(22b) 



which yields, in the Newtonian limit, an effective accel- 
eration of magnitude 
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(23) 



where if is a constant. Numerically, we obtain that the 



correction decays like r 



(for (i = 0.. 



-0.866 



(for 



(3 = 4.4) and r ~°- 762 (for (3 = 14). It is exactly Milgrom 
(oc 1/r) if (3 = 1. 

Lastly, we compute the rotation curves for our mod- 
ified dynamics, assuming a modified Hubble profile JTcj ] 
to find M{R). In figure [| we plot the rotation curves for 
K = 0.1 and several values of j3, as well as the curves for 
/3 = 4.4 for several values of K. The figure shows that, 
within the current limits on the Higgs mass, the modified 
dynamics allows for a broad family of curves, depending 
on the parameter K. Milgrom-type corrections are also 
allowed, although an exactly Milgrom 1/r dependence 
seems rather unlikely. 

Although more tests would clearly be needed, the ex- 
istence of a model such as the one presented in this Let- 
ter shows that one could imagine a situation, compatible 
with observations of the dynamics of galaxies and galaxy 
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FIG. 3. The coefficients ccr and ctg of the equations of state 
for the NO vortex, as functions of f3. 
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clusters, where dark matter would not be needed at all 
(or, taking a different point of view, would consist in our 
case of a fluid of vortices). This would have the advan- 
tage, over other dark matter models, that it is based on 
a particle theory (albeit a toy model here), giving it a 
stronger theoretical base. 
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FIG. 4. Sample rotation curves obtained from the acceler- 
ation (p^). On the top figure, we plot the curves for (bottom 
to top) (3 = 0.88, 1, 4.4 and 14. On the bottom figure, we plot 
the curves for f3 = 4.4 and (bottom to top) K = 0.1, 0.25, 0.5 
and 1. 
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